Gaussian quantum computation with oracle-decision problems 
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We study a simple-harmonic-oscillator quantum computer solving oracle decision problems. We 
show that such computers can perform better by using nonorthogonal Gaussian wave functions 
rather than orthogonal top-hat wave functions as input to the information encoding process. Using 
the Deutsch-Jozsa problem as an example, we demonstrate that Gaussian modulation with opti- 
mized width parameter results in a lower error rate than for the top-hat encoding. We conclude 
that Gaussian modulation can allow for an improved trade-off between encoding, processing and 
measurement of the information. 
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I. INTRODUCTION 

The process of computation involves taking an input 
and converting it into an output through the application 
of an algorithm [l|, 0] . There are two versions of classical 
computation: digital and analogue. Digital computers 
use finite-length bit strings as the input and output of the 
computer. Analogue computers typically use the electri- 
cal quantities of inductance and charge as analogues of 
mass and displacement so that continuously varying volt- 
ages may output, for example, the simulated trajectory 
of a spacecraft Q. 

Quantum computation also has two versions of infor- 
mation processing referred to as discrete-variable [J| and 
continuous-variable [5] quantum computation. In the 
discrete- variable version, algorithms employ a finite num- 
bers of qubits [I[ or qudits Q analogously to bits in the 
classical digital case. Continuous-variable quantum com- 
putation uses continuously-parameterized quantum sys- 
tems to process discrete information. 

Comparing the performance of a well-known problem 
and its algorithm in both discrete and continuous settings 
provides useful insight into how the different versions of 
quantum computation differ in their performance. The 
historically important Deutsch-Jozsa problem [3, H| has 
been used to show the different algorithm performance 
that can be achieved in the two settings. In the discrete 
variable setting, the Deutsch-Jozsa problem can solved 
exactly with a single query of a quantum oracle Q . In the 
continuous variable setting, quantum algorithms employ- 
ing orthogonal wave functions as the computational basis 
necessarily have single-query success probability less than 
one Q. 

Continuous-variable studies are often based on quan- 
tum optics because of the wide variety of tools that 
have been dev elop ed to process and measure optical field 
modes [f| [Iol - ll7l |. Continuous-variable quantum algo- 
rithms have been studied that use a single mode [9( and 
multiple modes fisij . In continuous- variable quantum in- 
formation procedures, the input state is typically a Gaus- 
sian wave function over the canonical position represen- 



tation with the physical system being a harmonic oscilla- 
tor (equivalently a single- mode optical field [HI)- How- 
ever the unbounded extent of these wave functions does 
not naturally fit with the finite length of the information 
strings being processed. 

To deal with the problem of finite-length informa- 
tion strings, the single-mode algorithm defined in Q has 
input states represented by orthogonal wave functions. 
Information is represented by finite-length bit strings 
z £ {0, 1}^ with N the number of bits. These strings are 
encoded into a region of the momentum domain p extend- 
ing from — P to P. A regular lattice of N discrete values 
of p are embedded in this domain such that the i th bit of 
z, is assigned to p as follows: p Zi := (I + 2i — N± 1)P/N . 
The lattice is thus {p Zi } with spacing 2P/N. 

An orthogonal basis of top-hat functions is formed 
from the p Zi , and the encoded momentum wave function 
is expressed as I^q 1 (-!)** \Pzt)- Tne ke ts \p Zi ) 
are phase-modulated by their corresponding bit values, 
and with this phase modulation, each of the possible 2 N 
strings is uniquely represented. Note that the constant 
wave function is the top-hat function extending from — P 
to P. 

There is a translational invariance between compu- 
tational basis states because each of the top-hat basis 
functions occupies an identically-sized region of momen- 
tum space and because the string z' is obtained from the 
string z by the translation z' = z © (z © z'). We can re- 
gard this finite basis as an infinite basis modulated by a 
top-hat function extending from —PtoP, which has the 
effect of truncating the allowed strings from an infinite 
domain to being from to N — 1. 

In this paper, we extend this approach by shaping the 
overall top-hat function with a Gaussian having width 
set by its standard deviation a. The Gaussian is trun- 
cated for \p\ > P, which has the effect that the com- 
putational basis states, although still orthogonal, are no 
longer translations of each other but are more compli- 
cated Gaussian-modulated basis states. At first glance, 
this latter feature would appear to provide a disadvan- 
tage, but using the mathematical properties of the error 
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function, we prove that the single-query success probabil- 
ity for the Gaussian case Pr^ > Pr^ is greater than the 
single-query success probability for the orthogonal case. 

This result is enabled by the extra degree of freedom 
manifest in the spread of the Gaussian wave function. 
Tuning the available parameters of encoding width, the 
spread of the Gaussian wave function and the width of 
the measurement window results in a more efficiently en- 
coded momentum wave function leading to the improved 
single-query success probability. Note that the use of 
the top-hat basis to encode information into a single 
harmonic oscillator is different than the approach used 
in where information is encoded into a collection of 
harmonic oscillators. 

Our paper is organized as follows. In Sec. UH we define 
oracle-decision problems and present the Deutsch-Jozsa 
0, HI, [13] in this context. We give an overview of the 
single-mode algorithm employing orthogonal states [9( 
that solves the Deutsch-Jozsa problem. We give a brief 
introduction into the coherent states of the harmonic os- 
cillator and define the single-mode algorithm in these 
terms. In Sec. IIII1 we prove that the single-query success 
probability claimed for the Gaussian model with trun- 
cated momentum is better than that achieved using or- 
thogonal states |9[. We conclude in Sec. IIVI 



the value of the function on the i bit string, taken in 
lexicographical order. 

For our analysis, we repeat the definition of an oracle 
decision problem given in [21| as follows. 

Definition 1. An oracle decision problem is specified by 
two non-empty, disjoint subsets A,Bc {0,1}^. Given 
a string z £ AU B = C , the oracle- decision problem is 
to determine whether z £ A or z € B with the fewest 
queries to the oracle possible. 

For completeness, we also repeat the definition of the 
Deutsch-Jozsa problem 0, [2(| in terms of Definition Q] 
given in [2l[ as follows. 

Problem 1. Given set the set of balanced strings A C 
{0, 1}^, where exactly N/2 elements take on the value 
and the set of constant strings B C {0, 1}^, where all 
N elements take on the same value everywhere, and a 
string z randomly selected with uniform distribution fi 
such that z C = A U B, the Deutsch-Jozsa Problem 
is to determine if z g A or z € B with the fewest oracle 
queries. 



B. Single Mode Continuous Variable Algorithm 
with Orthogonal States 



II. BACKGROUND 

An important aspect of quantum information process- 
ing is the ability to solve oracle decision problems with in- 
creased efficiency compared to classical information pro- 
cessing. In the case of oracle decision problems, efficiency 
is measured in terms of the number of oracle queries re- 
quired to solve the problem. Comparing the single-query 
success probability of an algorithm in both the discrete 
and continuous quantum settings provides useful insights 
into the advantages of one setting over the other. The 
Deutsch-Jozsa oracle-decision problem has been studied 
in both the discrete and continuous-variable quantum 
settings H, [t| HJ. Here we explore a single-mode 
continuous- variable quantum algorithm where the input 
state is a Gaussian wave function. 



A. Oracle Decision Problems and the 
Deutsch-Jozsa Problem 

The challenge of an oracle decision problem is to iden- 
tify which of two mutually disjoint sets contains a unique 
iV-bit string by making the fewest possible queries to an 
oracle. The oracle decision problem is typically couched 
in terms of a function / that maps N = 2™-bit strings to 
a single bit 



/: {0,1}" -+{0,1}. 



(1) 



Any Boolean function on n bits can also be represented 
by a string of N — 2" bits, in which the i th bit Z{ is 



In Fig. [T] we present the single-mode, continuous- 
variable quantum algorithm [9[ for the solution of oracle 
decision problems. The vertical lines on Fig. [T] repre- 
sent the states after the various steps of the algorithm 
using function notation rather than Dirac notation. In 
function notation, the Dirac ket \<p) is represented by the 
square-integrable function 



(x) = (x\ 



(2) 



where x in this case is the continuous position variable. 

The square-integrable condition means that orthogonal 
functions may be used to represent the wave functions. 
One possible set of orthogonal functions is the Fourier- 
transform pair realized by the sine/top- hat functions [9[. 
In this case, the sine function 



<f> {x) 



sin(Pa;) 



(3) 



is the input state, and its Fourier transform is the mo- 
mentum top-hat function 



i 



l, if p e [-P,P] 
0, i£p£[-P,P], 



(4) 



having finite extent of 2P in the momentum domain. One 
nice feature of the sinc/top-hat pair is that the finite 
extent of the top-hat distribution allows for finite length 
information to be encoded naturally. 

The encoded position sine function has unbounded ex- 
tent, and analysis of the optimum position measurement 
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itself. Therefore for notational simplicity, we chose to use 
the position representation of the vacuum, 



(x|0) 



(8) 



0o(x) 



0i (P) 



<fe(P) 



&00 



FIG. 1: Single-mode quantum circuit implementing the 
continuous- variable Deutsch-Jozsa algorithm [jj. 



window reveals an uncertainly relationship [i| between 
the measurement window S and the encoding length P 
expressed as 



PS 



7T 

2 ' 



(5) 



As a result, the single-mode continuous variable algo- 
rithm is necessarily probabilistic 9] and has single-query 
success probability 



Pr^ = 0.61. 



(6) 



Here we demonstrate that this single query success prob- 
ability may be improved upon by using a Gaussian wave 
function as algorithm input. 



C. Single Mode Continuous Variable Algorithm 
with Gaussian States 

The sine function employed as algorithm input in Q 
cannot be readily created in the laboratory. Here we are 
inspired by the ability to create and manipulate physical 
states of light in the laboratory using the tools of quan- 
tum optics. In particular, we employ coherent states, 
which may by represented by Gaussian wave functions, 
as the input states to our algorithm. 

The method of coherent states is well established and 
one feature is that the coherent states are overcom- 
plete [22|, (23[. High quality lasers generate light fields 
that are coherent [19j. The vacuum state is a displaced 
coherent state and as such has the same quantum noise 
properties. Coherent states are usually expressed as the 
ket | a) with a — XQ+ipo to reflect that it is a state that is 
shifted from the vacuum by the magnitude \a\. Similarly, 
the vacuum is usually expressed simply as the ket |0). 

In the position representation, the coherent state of 
laser light may be expressed as 

e -^[(x~x ) 2 -2ipox+ipox ] 

cf)(x; a) = {x\a) = -= , (7) 

where Xq = Po = corresponds to the vacuum state. 
Prom the perspective of our quantum algorithm, the dis- 
placed vacuum behaves no differently than the vacuum 



as the starting-point state for our algorithm. 

Quantum optics has many tools that allow for the ma- 
nipulation of light. Of interest in our algorithm is light 
squeezing, where quantum uncertainties are redistributed 
altering the shape of the distribution. The squeezing op- 
erator is given in as 



s(0 



(9) 



where a = x + ip is the annihilation operator and = 
x — ip is the creation operator. The quadratures x and 
p are regarded as the position and momentum of the 
harmonic oscillator, and the quantity Q is referred to as 
the squeezing parameter [19j . 

In Dirac notation, the squeezed vacuum state may be 
expressed as 



\<M = 5(C)|o). 



(10) 



In our analysis, we use the standard deviation a = 

to represent the effect of the squeezing operator on our 

function representation of a Gaussian state. 

We employ the squeezed vacuum as the input state to 
our algorithm, which we represent in function notation 
as 



<h(xiv) = (x\S(()\0) 



(11) 



The subscript zero identifies this state as the algorithm 
input state represented by the leftmost vertical line in 
Fig. (fTJ). We prove that the algorithm with this Gaussian 
input state has improved single-query success probability 
over the algorithm employing orthogonal states as input. 

Theorem 1. Using the single-mode quantum circuit 
given in Fig. [7] with the coherent state given by Eq. ill]) 
as algorithm input and employing sharp information 
cutoff, the single-query success probability Pr^ > Prj; 
is greater than the single-query success probability Pry- 
obtained using orthogonal states given by Eq. 



III. BOUNDING THE QUERY COMPLEXITY 
OF THE SINGLE MODE ALGORITHM WITH 
GAUSSIAN INPUT STATES 

The continuous-variable quantum algorithm using or- 
thogonal states solves Problem 1 with exponentially 
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small error probability 1/(2™) in a linear number of 
queries 6(m) @. This query complexity is dependent 
on the single-query success probability Pr ^ , which is a 
measure of the maximum achievable separation between 
the probability that the encoded string is a balanced 
string versus a constant string. Here, where the input 
is a Gaussian state, we demonstrate that the key param- 
eters affecting this separation are the encoding width, the 
spread of the Gaussian wave function and the width of 
the measurement window. We vary these parameters and 
discover their optimum values in our proof of Theorem 1 . 



A. Encoding Information into Gaussian States 

With reference to Fig. HI the first step of the algorithm 
is to take the Fourier transform Q of the input state 
<t>o(x\o) giving 



12 2 



(12) 



The next step has the oracle U / modulate the momen- 
tum Gaussian with the pulse train that represents the 
encoding of the iV-bit string z. 

The modulated momentum wave function is 



hip) = C v > ;o,p) = ftfWfc p)Mp ; *), (13) 

where we have labelled the state with all relevant pa- 
rameters. Descriptions of the elements of this equation 
follow. The modulating square-wave encoded with the 
TV-bit string z is 



N-l 



/W(p;P) = ^(-l)-nW( P ;P), 



(14) 



The encoded position state is thus expressed as 



(x) = 4,W(x;a,P) 



Mi N) (x;a,P). (18) 



The effect of the encoded information is completely cap- 
tured in the position modulating term 



N 



MW(x;a,P)=^(-l) 
4=1 

— erf 



erf 



•djcr 2 + ix 
V2a 



where 



-dj-ia 2 + ix 

2j- A T 
N 



(19) 



(20) 



The final step is the measurement step. 

We follow the same approach taken in Q and calculate 
the probability of detecting a particular wave function in 
the interval ±<5 as 



<t>W(x;a,P) 



(x;a,P) 



dx. 



(21) 



Since the wave function may be encoded with a constant 
string or a balanced string, we need to determine the op- 
timal value of S that maximizes our ability to distinguish 
between these cases. Our approach is to determine which 
balanced functions dominate all other balanced functions 
in the measurement window. 

We begin by defining three pairs of iV-bit strings: the 
antisymmetric balanced (AB) strings, the symmetric bal- 
anced (SB) strings and the constant (C) strings as 



where the definition of the momentum bins given in Q 
is repeated here as 



1 0, otherwise. 



N 



2^4=*), -(1-2^)] 
(15) 



Note that the modulating function has the effect of chop- 
ping off the tails of the momentum Gaussian outside ±P 
thus truncating the Hilbert space. 

The normalization factor, rj, of the chopped distribu- 
tion is calculated as 



& n) (p;<t,p) 



dp = crf(Pcr) 



(16) 



where the error function is erf(ui) = f™ e 1 dt, and 

r/ = l/y/ed(Pa). (17) 

The penultimate step is to take the inverse Fourier trans- 
form of this encoded momentum state. 



AB G < 0-.-01- • ■ 1, 1- ■ • 10 - • - 0. 

N/2 N/2 N/2 N/2 



(22) 



SB G < .0 ■ ■ ■ 01 ■ ■ ■ 1.0 ■ ■ ■ 0, 1 ■ ■ ■ 1.0 ■ • ■ 01 ■ • ■ 1 > , (23) 

N/4 N/2 N/4 N/4 N/2 N/4 



C6{0-0,1-1 

TV N 



(24) 



Note that the constant stings have zero bit transitions, 
the antisymmetric balanced strings have one bit transi- 
tion, and the symmetric balanced strings have two tran- 
sitions. All other balanced strings have two or greater 
transitions. 

It is insightful to analyze the modulating term given 
by Eq. (fl~9| for x — 0, which we express as 



N 



Mf)(0;a,P)=^(-l) ! 

4=1 



erf 



(25) 
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We use the anti-symmetric property of the error function 
erf(a) = —erf (—a), and the property that erf(O) = 0. Wc 
also use the facts that i9jv = P, = and dj = 

— $N-j for j = 0, 1, . . . , N in determining the following 
results. 

For the constant case, all the terms cancel except the 
first and last, and we have 



M^(0;a,P) = ± 



erf 



$NCr 
V2 



MS- 



(26) 



For the antisymmetric balanced case, all the terms cancel 
and we have 



M% ) AB (P;<r,P) = ± 



erf 



V2 



erf 



\V2j 



(27) 



For the symmetric case we have 




FIG. 2: Plots of the magnitude of the position modulation 

function |m z (8) (as; a, 1)1 given by Eq. {TSJ for (a) a = 0.4, (b) 

a = 0.6, (c) a = 0.8, and (d) a = 1.0. For |x| < x c , the 
symmetric balanced function (bold) dominates all other bal- 
anced functions, and for |as| > x c , the antisymmetric balanced 
function (bold) dominates all other balanced functions. 



M^l B (0;a,P) = ± 



(28) 



Here the sum is non-zero except for in the limiting case 
where Per — > 0. 

It is apparent from the results of Eqs. (|27p and (|28|) 
that there are different classes of balanced functions. 
Some balanced functions are only non-zero at x = in 
the limit as a goes to zero, and some balanced functions 
are zero at x — for all values of er. We need to de- 
termine which balanced functions are the two functions 
that dominate the measurement region. 

In Appendix A, we prove that the magnitude of 



M { z N \x;a,l) 



given 



the position modulation function 
by Eq. (|T8|) and subject to the balanced condition 
2fcLi( — l) Zi = 0, is maximized by either the antisym- 
metric balanced function given by Eq. (|22[) or the sym- 
metric balanced function given by Eq. (|23|) . For N = 8, 
the situation is presented in Fig. [31 where it be seen that 
the actual dominating function is dependent on the value 

of (J. 

In Fig. [3J the crossover point is drawn and is approxi- 
mated in Appendix [A] as 



Pa 2 



(4-PV 2 )' 



(29) 



For | a: | < x c , the symmetric balanced function (shown in 
bold) dominates, and for |x| > x c , the antisymmetric bal- 
anced function (shown in bold) dominates. All remaining 
balanced functions, of which there are a total of (?) = 70 
are depicted as light gray lines in Fig. [H We use these 
results to complete the proof of Theorem 1. 



B. Proof of Theorem 1 

We need to maximize the separation between detecting 
a balanced string and a constant string. To this end, we 
define the following quantities 



(<?,P,5) 



and 



Pr 



Pr 



9 zGC 



l(N) 



Pr 



Pr 



; z£AB 



-zssb 



(30) 



(31) 



where for brevity, we have suppressed the arguments in 
Eq. (|2ip . The single-query success probability is defined 
in these terms as 



Pr^ = min [A AB (a, P, 6) , A SB (a, P, 5)} 



(32) 



This expression assumes that either the antisymmetric or 
the symmetric balanced strings dominate all other bal- 
anced strings in the region ±5 as presented in the pre- 
vious subsection and proved in Appendix A. We seek to 
determine the values of 8 and a that maximize the sepa- 
ration between these two probabilities. 

With A AB (S,a,P) defined in Eq. (JSOJ), we 
set A' AB (S : cr,P) = 

A' AB (S,a,P) = 



JjA A b(<5, cr, P), which gives us 



20Faerf(Pcr) 
■Pa 2 - 



2 erf 



2 erf 



erf 



V2a 
Pa 2 + 



f Pa 2 -\8 
V V2ct 

(4i 



erf 



/Per 2 +iS 



erf 



V\/2cr 
S 



+ 2 erf 



( ' lS 



V2c 



-Hi 



V2o 



(33) 
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It suffices to set A' AB (5, cr, P) — to maximize the sepa- 
ration. Before doing so, we elect to simplify Eq. (j3"3")l by 
'normalizing' the standard deviation a and the measure- 
ment 'length' 5 with respect to the encoding 'length' P. 

We assume that the uncertainty relation [9j remains 
true up to a constant. We express this as 



PS = S. 



(34) 



This assumption and analysis of the error function argu- 
ments of Eq. p3[) result in a similar uncertainty relation- 
ship between P and cr, which we express as 



(35) 



Pa 



Making the substitutions given by Eq. (|34[) and Eq. ([35]) 
into Eq. (|33p and setting it to results in the following 
expression 



A^ B (S,a)=0 



erf 



2 erf 



2 erf 



i<5 



V2a 



erf 



iS 



V2a 



iS 



y/2a 

T 2 +iS 



erf 



( [S 



V2d 



V\/2> 

-(£) 



2 erf 



(36) 



Note that the variables cr and 8 are, in some sense the 
'normalized' Gaussian standard deviation a and the mea- 
surement width 5, 'scaled' by the momentum 'length' P. 

Eq. (|36[) is dependent on the two variables, 6 and cr and 
is thus insufficient to find the global optimum values of cr 
and <5. We obtain the needed constraint from the similar 
equation derived from the symmetric balanced function. 
Following the same steps we did in Eq. (|33)l and Eq. (|36|) , 
we obtain the following expression 



Kb (MHO 
= erf 



-erf 



erf 



2y/2a 



\2V2aJ 



a 2 + iS 
2^2a 



erf 



erf 



erf 



2^2a 
<r 2 -i~8 

2 +i5 



V25 



(37) 



We solve Eqs. ([31)]) and (|3"T1) simultaneously to establish 
the optimum values of the measurement lengths i5ab and 
(5sb in terms of the normalized standard deviation a. 

In Fig.[31 we plot the distributions \4>i N \x; cr, P)\ 2 for 
z G {AB, SB, C} for several values of a. We also plot ver- 
tical lines corresponding to the values of 5 corresponding 
to A AB = and A SB = 0. Note that there are values of 
the normalized parameters a and 5 where A AB (<5, cr) = 0, 
and A SB (<5, ct) = simultaneously. This situation occurs 
where 6 « 2.30 and a w 2.11 and is depicted in Fig.^b). 
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(d) 
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-3 -1.5 



1.5 



FIG. 3: Plots of ^(icjer,!) for z e {AB, SB, C} (solid, 

dotted, dashed respectively) for: (a) a — 1.67, (b) a = 2.11, 
(c) a — 2.5, and (d) a — 3.0. The solid vertical lines in all 
four sub-plots correspond to A' AB (5, a) = 0, and the dashed 
vertical lines correspond to Ag B (<5, a) = 0. Note that a) cor- 
responds to the values that optimize the single-query success 
probability. 



However, these values do not optimize the success 
probability since 



A AB (2.30,2.11) 0.68, and 
A SB (2.30,2.11) 0.54. 



(38) 



Lack of optimality is manifest in the lower of the two 
above values, which is less than single-query success prob- 
ability for the orthogonal case Pr^ pa 0.61. 

Increasing the value of a further serves to in- 
crease A AB (5, a) and decrease A SB (S, a) , which worsens 
the success probability. Reducing the value of a brings 
them together. The quantity A AB (<5, cr) thus takes on its 
maximum value when 



subject to the constraint 

A AB (6, a) = A SB (6,a) . 

2.01 and a 



This occurs at a value of <5ab ss 
Optimality is manifest since 

A AB (6, a) = A SB (5, a) 



0.68. 



The optimal situation is depicted in Fig.[3][a). 
For P = 1, we express the optimal parameters as 



and 



2.01, 



1.67. 



(39) 

(40) 
1.67. 

(41) 

(42) 
(43) 
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FIG. 4: Comparison between the orthogonal (solid) and the 
Gaussian (dashed) wave functions. The amplitude of the en- 
coded momentum functions for the constant case is depicted 
in (a). The respective algorithm success probabilities are de- 
picted by the shaded regions of (b), where light gray corre- 
spond to the orthogonal case and dark gray corresponds to 
the Gaussian case. Medium grey corresponds to the region 
of overlap. Note that the encoded momentum Gaussian is on 
average narrower than the orthogonal pulse, which results in 
respectively broader position wave functions. 

For these values, the single query success probability of 
the Gaussian model with sharp information cut-off model 
is 

PrJ, « 0.68. (44) 

This upper bound is approximately 10% greater than 
that shown for the model employing orthogonal states 
0, where Pr^ w 0.61 and 6 X = § . 

At first glance, the increase in single-query success 
probability of the Gaussian case over the orthogo- 
nal case appears somewhat surprising. The Gaussian 
wave functions are coherent states and therefore non- 
orthogonal [HI . Intuitively, the orthogonal states should 
be optimal especially given that the finite extent of the 
momentum wave functions provides a natural fit for en- 
coding finite infirmation. 

Upon closer inspection however, we see that the im- 
provement results from the ability to 'tune' the Gaus- 
sian spread, represented by er, to match the encoding 
length P. No such 'tuning' is possible with the finite 
states. We depict this in Fig. \M,&) for the constant case 
with P = 1 and optimal cr" = 1.67. We see that the en- 
coded momentum Gaussian wave function is on average 
narrower than the orthogonal pulse wave function. Since 
the momentum and position wave functions are Fourier 
transform pairs, narrowing of one results in broadening 
of the other. 

The subsequent broadening of the encoded Gaussian 
wave functions results in a wider optimal measurement 
window <$" > 6^. This leads to a greater single-query suc- 
cess probability and is represented by the shaded regions 
in Fig. H^b). The larger dark gray region corresponds to 
the single-query success probability offered by the Gaus- 
sian wave functions. We thus conclude that the increased 
success probability is achieved through the extra degree 
of freedom afforded by a. For P = 1, this requires that 
the input state be squeezed to a sa 1.67. 



IV. CONCLUSIONS 

We have shown that a simple-harmonic-oscillator 
quantum computer solving oracle decision problems per- 
forms better using non-orthogonal Gaussian wave func- 
tions as the algorithm input rather than the orthogonal 
top-hat wave functions. We have also shown that the lim- 
iting case of the Gaussian model for a — > and non-zero 
P corresponds to the model employing orthogonal states. 
In both cases, the computational bases are orthogonal, 
and encoding takes place in the momentum domain and 
information processing and measurement take place in 
the dual position domain. Also in both cases, the single- 
query success probability is dependent on the maximum 
separation between the position wave function encoded 
with the constant string and the position wave function 
encoded with the worst-case balanced string, which is the 
antisymmetric balanced string. 

In the orthogonal case, ./V-bit strings are uniquely en- 
coded into the computational basis formed by the top-hat 
functions, and the overall-width of the encoded string is 
set by the encoding length P. In the dual position do- 
main, the encoded string is represented by a sum of N 
equi- angularly spaced, equi-length phasors multiplied by 
a sine function. The rate at which the constant sine func- 
tion falls off its peak and the rate that the antisymmetric 
balanced sine function rises from its minimum sets the 
size of the optimum position domain measurement win- 
dow. Thus the optimum position domain measurement is 
set by sharpness of the sine function, which is dependent 
on the encoding length only. 

In the Gaussian case, TV-bit strings are uniquely en- 
coded into the computational basis formed by more com- 
plicated Gaussian-modulated basis states. The overall- 
width of the encoded string is again set by the encoding 
length P, but it is also shaped by the Gaussian spread a. 
In the dual position domain, the encoded string is rep- 
resented by a sum of non equi-angularly spaced and non 
equi-length phasors multiplied by a Gaussian function. 
The rate at which the constant encoded function falls off 
its peak and the rate that the antisymmetric balanced 
function rises from its minimum is governed by both the 
Gaussian spread a and the encoding length P. More im- 
portantly, the rate set by the optimal values of a and 
P is more gradual than that achievable in the orthogo- 
nal case allowing for greater separation between the two 
probabilities. 

We thus conclude that the Gaussian allows for an im- 
proved trade-off between encoding, processing and mea- 
suring of the information. Encoding takes place in the 
momentum domain, and the Gaussian takes better ad- 
vantage of the space available to encode the informa- 
tion. Correspondingly, information processing and mea- 
surement take place in the dual position domain. The 
Gaussian-encoded position wave function enables a wider 
measurement window, which means more of the encoded 
information is available for distinguishing between a wave 
function encoded with a constant string and a wave func- 
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tion encoded with the worst-case balanced string. 
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Appendix A: Proofs of dominance of symmetric and 
antisymmetric balanced functions 

In this appendix we prove that the symmetric and the 
antisymmetric balanced functions maximize the magni- 



tude of 



M { z N) (x;a,r 



given by Eq. (JT5J) subject to the 



balanced condition ^2f =1 {— l) Zl = in three Lemmas. 
Lemma 1. For the region \x\ < tt with a — and 



subject to the balanced condition 53j=i( — I) 2 
max (j)i N \x; a, P) occurs for z G {AB}. 



0. 



Proof. We prove this Lemma by showing that, in the lim- 
iting case where a — > 0, the encoded position wave func- 
tion given in Eq. (fT8"|) becomes the same as the encoded 
orthogonal wave function analyzed in [9[ . In that case, it 
is proved that the antisymmetric balanced function dom- 
inates all other balanced wave functions in the region of 
interest. 

We begin by defining the quantity 



A k (x,a) = erf | 



2Pk 2 

N 



IX 



V2o 



/ ap(fc-i) a 2_^ 
-erf - N 



V2a 



(Al) 



for k = I,..., N/2. We use this term here and in later 
Lemmas. For ease of understanding the antisymmetric 
features of this term, we have elected to change the count- 
ing variable in the term §j given in Eq. (|20l) from j to k, 
where j = k + N/2. We express the k th term of the en- 
coded position wave function given by Eq. (|18j) in terms 
of this quantity as 



^ z k \x-a,P)=± 



rj e 



A k (x,a) 



(A2) 



where the ± represents the effect of the bit (— l) Zfe . We 
represent this quantity as the phasor 

4 k) {x; a, P) = ±R k (x; a, P) exp [i<p k {x;a,P)], (A3) 

to align with the description of the orthogonal case. 



The phasor magnitude is expressed 

x 2 

r\ e 



Rk(x;a, P) 



2^J5 



A, 



and the argument is 



tpk{x;a, P) = arctan 



/ Im [A k ] 
\Re [A k ] 



(A4) 



(A5) 



where we have suppressed the arguments of A k for the 
sake of brevity. 

The quantities R k (x;a) and ifk(x;o~) are too opaque 
to understand limiting behaviour, so we use Taylor series 
analysis to gain insight. The Taylor series representation 
of angle ipt(x; a) given by Eq. (|A5|) is expressed 



{2k-l)x (2k-l)xa 2 2(2fc-l)xcr 4 



N 37V 3 
+ O (xa 6 ) + O (xV) 

where for a = 0, we have 



45A 5 



(2k-l)x 

N : 



(A6) 



(A7) 



which presents an equiangular separation between subse- 
quent phasors. 

Similarly the Taylor series for magnitude given by 
Eq. (|A4[) is expressed 



R k (x;a) 



1 



(N 2 - 12(/c - ~ 4 ) °- 2 
N\/^ " 6AV^ + 

(-57V 2 + 60(fc- l)k + 24)a 2 x 2 n . , 4N 

\-(J [x a ) . 

(A8) 



18OA 5 0F 



For cr = 0, this gives 
VP 

Rk(x;0) 



p5/2 x 2 



2 m 



6A 3 0F 
sin(Px/A) 
V Pirx 



(-iryp(^) 

A0F(2m+ 1)! 

(A9) 



where the last step assumes the limit m — > oo. 

Combining the results of Eqs. (fAljj) . (|a7|) , and (K9\i 

gives 



e)(,;0,P) = ^fe^£(-l)^e i ( 
•'Pttx 



N-(2j-l) 



)Px 

(A10) 



which is the encoded orthogonal position wave function 
given m @. The proof given therein suffices to prove 
Lemma 1. □ 

Lemma 2. For x — and a > and subject to the bal- 



anced constraint X)j=i(~l)^ J = 0; max 
curs for z £ {SB}. 



S?°(0; a) 



oc- 
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Proof. We exploit the structure of the quantity given in 
Eq. (|A1|) with x — expressed as 



erf 



(All) 



Using the shorthand = A k (0,a), we express a set of 
TV terms in the following convenient form 



{A,,.. 



,A 2 ,A 1 ,A 1 ,A 2 , 



,Ai 



(A12) 



Note that the A k are real numbers. 



for k = 1,2, 

We now show that Afc > v4fe +1 . We express the differ- 
ence between these terms as 



A k - A 



k+l 



2 erf 



V2kc 
N 



erf 



'V2(k-l)a 



N 



erf^ 



Da 



N 



(A13) 



Showing that Eq. (|A13|) is positive for all k is equivalent 
to showing that 



2 erf(a) - erf(a + b) - erf(a - b) > 



(A14) 



for a, & € ffi. and a, 6 > 0. 

Over the domain (0,oo), the error function is strictly 
monotonically increasing with strictly monotonically de- 
creasing slope g^erf(x) = 2e~ x "/y / 7r. This means that 
successive increments Sx result in decreasing Sy = erf(<5x) 
increments. This may be expressed as 

erf(a) - erf(a - b) > erf(a + b) - erf(a), (A15) 

and thus 

2erf(a) - erf(a - b) - erf(a + b) > 0, (A16) 

which establishes that Ak > A k+ i. 

The strategy required to maximize the sum of the N 
terms of the set (|A12[) subject to the balanced constraint 
is now clear. Since A\ > A 2 > ■ ■ ■ > An_ , the maxi- 

2 

mal term must contain as many of the larger terms as 
possible. This maximal sum is thus expressed 



±2 (a x + A 2 + • • ■ + A, 



Ai 



A 



£ +2 



A 
(A17) 



This expression manifests the symmetric balanced (SB) 
function definition given in Eq. (|23[) thus proving the 
Lemma. □ 

Lemma 3. For \x\ > and a > and subject to the 



balanced condition J2j=i(~ = 0, max <j>z (x;a) 
curs for either z € {SB} or for z € {AB}. 



oc- 



Proof. We modify the set of elements Ak to include the 
imaginary components resulting from x > as 

|v4«_(a;, a), . . . , A\{x, a), A±(x, a), . . . , An_(x, cr) } . 

(A18) 

We now exploit the antisymmetric property of this set. 
The fact that erf(u>*) = erf*(ui) allows us to use the 
notation 



and 



A k (x,a) = a k + i/3 k , 



A* k (x, a) = a k - i/? k 



(A19) 



(A20) 



to capture the overall of effect of the error function having 
complex arguments. 

The strategy to maximize the sum of the elements in 
the set given expression (|A18[) subject to the balanced 
constraint is clear. The sum must be either purely real or 
purely imaginary. A complex sum reduces these achiev- 
able maximums in two ways. It causes elements to be 
subtracted, and it results in a vector sum rather than a 
liner sum. 

The maximum real sum subject to the balanced con- 
straint is 



k=N/4 



k=N/2 



J2 A k (x,a)+A* k (x,a)- ^ A k (x, a) + A* k (x, a) 



k=l 

/ k=N/4 

= 2 i ^ a k (x,cr) 

k=l 




(A21) 



which is achieved for the symmetric balanced function 
demonstrated in Lemma 2. The maximum imaginary 
sum subject to the balanced constraint is 



k=N/2 

E 

k=l 



k=N/2 

A k (x,a)- A *k(%,<r) 
fe=i 

k=N/2 

= 2i ^(x,^, 

k=l 



(A22) 



which is achieved for the antisymmetric balanced func- 
tion. 

As x increases from zero, the imaginary component 
of the error function increases accordingly. For small x, 
the real part still dominates and the symmetric balanced 
function is the balanced function with the greatest mag- 
nitude. However, there is a point where the antisymmet- 
ric balanced function takes over the dominate role. We 
determine the value of this crossover point, x c , in terms 
of cr and P in the following. 

The N — 4 case is the simplest case which demon- 
strates the crossover. For this case the set is 



{A* 2 (x, a), A\ [x, cr), A x (x, a),A 2 {x, a)} . 



(A23) 
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The antisymmetric balanced sum is 



for which the lowest-order Taylor approximation is 



(-A*{x,a) - A\{x,a) + A\{x,o) + A 2 (x,a)) 

= a-i + a 2 + i(/3i + (3 2 ) - «i - a 2 + + /3 2 ) 

= i2(ft+^ 2 ), (A24) 

and the symmetric balanced sum is 

{-A*(x,a) + A\{x,a) + A^x, a) - A 2 (x, a)) 

= 2(ai - a 2 ). (A25) 

The switch over thus occurs when 

Wi+/3 2 ) = (a 1 -a 2 ), (A26) 



Pa 2 



(4-P 2 a 2 )' 



(A27) 



For N = 8, this crossover point from symmetric to anti- 
symmetric dominance is plotted in Fig. [21 □ 
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